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ABSTRACT 

In this work, we investigate if gravitational microlensing can magnify the polarization signal 
of a stellar spot and make it be observable. A stellar spot on a source star of microlensing 
makes polarization signal through two channels of Zeeman effect and breaking circular sym¬ 
metry of the source surface brightness due to its temperature contrast. We first explore the 
characteristics of perturbations in polarimetric microlensing during caustic-crossing of a bi¬ 
nary lensing as follows: (a) The cooler spots over the Galactic bulge sources have the smaller 
contributions in the total flux, although they have stronger magnetic fields, (b) The maximum 
deviation in the polarimetry curve due to the spot happens when the spot is located near the 
source edge and the source spot is first entering the caustic whereas the maximum photometric 
deviation occurs for the spots located at the source center, (c) There is a (partial) degeneracy 
for indicating spot’s size, its temperature contrast and its magnetic induction from the de¬ 
viations in light or polarimetric curves, (d) If the time when the photometric deviation due 
to spot becomes zero (between positive and negative deviations) is inferred from microlens¬ 
ing light curves, we can indicate the magnification factor of the spot, characterizing the spot 
properties except its temperature contrast. The stellar spots alter the polarization degree as 
well as strongly change its orientation which gives some information about the spot position. 
Although, the photometry observations are more efficient in detecting stellar spots than the 
polarimetry ones, but polarimetry observations can specify the magnetic field of the source 
spots. 


1 INTRODUCTION 

Gravitational microlensing is an astronomical phenomenon in 
which the light of a background star is magnified due to passing 
through the gravitational jield of a foreground object by producing 
distorted images l iEinstein 1936h . This phenomenon has many ap¬ 
plications in detecting extra-solar planets, measuring the mass of 
isolated stars, studying the stellar atmospheres, etc. (see, e.g. Mao 
2012,Gaudi 2012). 

One of features of gravitational micr olensin g is th at i t causes 
a net polarization for s ource stars llSchne ider & Wagoner 19871 : 
[Simmons et al. 1995a,bi : iBogdanov et al. 199(f) . There is a local 
polarization over th e sta r surface due to the electron scattering 
in its atmosphere l lChandras ekhar 196(J). Since, these local polar¬ 
izations have symmetric orientations with respect to the source 
center, the net polarization of a distant star is zero. Some effects 
can break this symmetry and create a net polarization for distant 
stars such as circum-ste ll ar disk of dust grai ns around the central 
star dDrissen et al. 19891 : lAkitava et al, 20091) . stellar spots, mag¬ 
netic field and lensing effect. During a microlensing event, the 
circular symmetry of source star is broken and result is a net po¬ 
larization for the source star. Measuring polarization during mi¬ 
crolensing events helps us to evaluate the finite source effect, the 
Einstein radius and the limb-darkening parameters of source star 
jYoshida 20061:1Agol 199&lSchneider & Wagoner 19871) . 

In single microlensing events, the polarimetric signal reaches 


to its maximum amount, i.e. about one percent, when the impact 
parameter of the source center i.e. u cm reaches to ~ 0.96p* where 
p+ is t h e projec ted source radius, normalized to the Einstein radius 
dSchneider & Wagoner 19871) . But, the probability of happening a 
microlensing event with an impact parameter as smaller as u C m < 
pi, is low. Whereas, in binary microlensing events whenever the 
source star crosses the caustic curve the polarization signal as high 
as one percent can be obtained dAgol 19961) . Indeed, the point-like 
caustic in a single lens converts to some closed curves in binary 
lensing which is extended about the angular Einstein radius (see 
e.g. Schneider and Weiss 1986). 

Consequently, whenever the source crosses the caustic curves, 
due to the high gradient of magnification and polarization, the 
anomalies over the source surface can be resolved through pho¬ 
tometric or polarimetric observations. One of the second-order per¬ 
turbations over the source surface is the stellar spot. Indeed, a mag¬ 
netic field over the source surface generates convective motions 
which make a temperature contrast with the source surface tem¬ 
perature, so-called as a stellar spot (e.g. Berdyugina 2005). A stel¬ 
lar spot produces a polarization signal through two channels (i) 
Zeeman effect due to its magnetic field and (ii) circular symme¬ 
try breaking of the source surface brightness due to the tempera¬ 
ture contrast with the host star. Comparing these two effects, the 
first one has a dominant contribution in polarization signal. A stel¬ 
lar spot over the source surface disturbs the light and polarimetric 
curves of microlensing events. Detecting and characterizing stel- 
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lar spots through their effects on the photometric light curves dur¬ 
ing microlenising events due to a point-like lens and binary lenses 
were extensively studied in (Heyrovsky & Sasselov (2000),Han et 
al. (2000), Chang & Han (2002), Hendry et al. (2002), Hwang & 
Han (2010)). Here, we investigate detecting and characterizing stel¬ 
lar spots over the source surface through their effects on the polari- 
metric curves in microlensing events. 

Although, most spot-induced polarization signals from the 
Galactic bulge sources are too weak to be detected by themselves, 
the lensing effect can magnify these signals and makes them be ob¬ 
servable. Studying spot-induced perturbations in polarimetric mi¬ 
crolensing due to a point-like lens was first investigated by Saja¬ 
dian and Rahvar (2014). They also noticed that there is an orthogo¬ 
nal relation between the polarization orientation and the astrometric 
centroid shift vector in microlensing events due to a point-like lens 
and binary lenses except near the fold of caustic and investigated 
the advantages of this orthogonality for detecting the anomalies of 
the source surface. Here, we study detecting and characterizing the 
stellar spot through polarimetric microlensing in caustic-crossing. 
We first explore the characteristics of spot-induced perturbations on 
the polarization degree and its orientation in binary microlensing 
events. Then, by performing a Monte Carlo simulation from binary 
microlensing events with spotted sources during caustic-crossing 
we estimate the probability of detecting a stellar spot through po¬ 
larimetric microlensing during caustic-crossing and the number of 
events. 

In the following section, we study the stellar spot effects on 
(i) magnification factor, (ii) polarization degree and (iii) polariza¬ 
tion orientations. In section 0 the detectability of the polarimetric 
and photometric deviations due to the stellar spots in microlensing 
events during caustic-crossing is studied by performing a Monte 
Carlo simulation. We explain the results in the last section. 


2 SPOT-INDUCED PERTURBATION IN 
POLARIMETRIC MICROLENSING 


In order to describe the polarized light we use the Stokes param¬ 
eters of Si, Sq, Su and Sv- These parameters describe the total 
intensity, two components of linear polarization and circular p olar- 
ization over the source surface, respectively iTinbergen 19961 ). The 
polarization degree ( P ) and the polarization an gles (j> v and_# p as 
functions of Stokes parameters are given by Jchandrasekhar 19600 : 



By considering the linear polarization of light scattering in a stel¬ 
lar atmosphere, the circular polarization is zero Sv = 0. In mi¬ 
crolensing events, the Stokes parameters during the magnification 
are given by: 
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Figure 1. 1 (//) the total intensity (solid red line), I— the polarized intensity 
(green dashed line) normalized to the source intensity at the disk center, i.e. 
/+, and the polarization degree (blue dotted line) versus the distance from 
the source center normalized to the source radius i.e. p. 



where p* is the projected source radius which is normalized to 
the Einstein radius, p is the distance from the centre of the stel¬ 
lar disk to each element over the source surface which is nor¬ 
malized to the source radius, p = y/l — p 2 , (j> is the azimuthal 
angle between the lens-source connection line and the line from 
the center of coordinate to each element over the source surface, 
u = («cm + p 2 pi + 2pp*w cm cos 0) 1 ' 2 is the distance of each 
projected element over the source surface with respect to the lens 
pOSltlOIl, Ucm is the impact parameter of the source center and mag¬ 
nification factor for simple microlensing is A(u ) = u , +2 . The 

uyj u 2 + 4 

amounts of I (p) and /_ (p), as the total and polarized light inten¬ 
sities, by assuming the electron scattering in spherically isotropic 
scattering atmosphere of an early-type star on the emitted light, 
were evaluated by Chandrasekhar (1960) as follows: 

I{p) = /*(l-ci(l -p)), 

I-(p) = I*c 2 (l-p), (3) 

where /* is the source intensity at the disk center, ci = 0.64 and 
C 2 = 0.032 dSchneider & Wagoner 19871) . Here, we focus on the 
giant sources with finite size effect for which the limb-darkening 
coefficients are d ifferent a bit. For these sources, we set ci = 0.7 
dZub et al. 20091) and C 2 = 0.04. 

Figure 0 represents I(p) (red solid line), I-{p) (green 
dashed line), normalized to the intensity of the source at the disk 
center, versus the distance from the source center normalized to the 
source radius i.e. p. According to this figure, considering the limb- 
darkening effect the total intensity I (p) decreases a bit whereas the 
polarized intensity /_ (p) increases from the center to the limb. If 
there is no lensing effect i.e. A(u) = 1, the local polarization de¬ 
gree is equal to P = I-(p)/I(p). Here, over the star surface the 
polarization degrees enhance from the center to limb as shown in 
Figure 0 with blue dotted line. Also, the polarization vector in 
each point of the source surface is tangent to a circle whose cen¬ 
ter coincides to the source center and passes from that point. The 
overall polarization signal of an ideal source is zero, because of the 
symmetric orientations of local polarizations with respect to the 
source center. 

Now, let us consider a stellar spot on the source surface. To 
specify the spot in our simulation, we use four parameters for spots 
as: (i) the size of spot, (ii) the location of spot on the source star, (iii) 
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its magnetic field and (iv) the temperature contrast of the spot with 
respect to the host star. These parameters are not completely in¬ 
dependent and there is a correlation relation between the magnetic 
field of t he spot and its te mperature contrast with respect to the 
host star feerdvugina 20051) . Also, we consider the following limi¬ 
tations for modeling the spot polarization: (i) There is one circular 
spot on the source surface, (ii) the magnetic field over the spot has 
a constant amount and its direction is radial, (iii) the magnetic field 
is high enough to ignore the Hanle effect (see e.g. Stenflo (2013)), 

(iii) we take the Minle-Eddington model for the stellar atmosphere 
in which the line source function is linear versus the optical depth, 

(iv) the ratio of the line and continuous absorption coefficients is 
constant, (v) the intensity over the spot is constant. Also we as¬ 
sume a black body radiation for the spot and its host star to indicate 
their intensities and finally (vi) we estimate the broad band Stokes 
intensities by integrating over the wavelength pass band, taking into 
account the Stokes intensities in a spectral line. 

The radius of the spot is r 3 and the angular radius of spot in 
coordinate located at the centre of source star is given by do = 
sm~ x {r s / R*). We preform a coordinate transformation to obtain 
the position of spot on the lens coordinate set two axes of which 
describe the sky plane. In this regard, we use the two rotation angles 
of 6 S around y-axis and <p 3 around 2 -axis. The more details about 
modeling of the spot in our simulation can be found in Sajadian and 
Rahvar (2014). 

To calculate the modified Stokes parameters, we integrate over 
the source surface from the Stokes intensities. Since, the Stokes 
intensities over a given spot are different from those over the other 
points of the source surface, we take into account the combination 
of the terms from spot and source star, using a step function. The 
overall Stokes parameters S'q, S'u, S' v and Sj are given by: 
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where 0(p, (j>) is that step function which is equal to one within 
the spot area and zero for the other points. Iq, 3 , Iu,s, Iv,s and Jr, s 
are the Stokes intensities over the spot which are measured accord¬ 
ing to the spot position, it magnetic field and its intensity. Indeed, 
the magnetic field through the Zeeman effect on the spectral lines 
causes a net polarization for the spot which is mostly circular (see 
e.g. Illing et al. 1975). The more details about calculating these 
Stokes intensities by considering the mentioned limitations over the 
source and its atmosphere are brought in Appendix section. 

The spot-induced perturbations on the magnification factor, 
polarization degree and polarization angles are explained in the fol¬ 
lowing subsections respectively. 


2.1 Perturbation on the magnification factor 

The spot-induced perturbation on microlensing light curve due to 
a point-like lens was first investigated by Heyrovsky and Sasselov 
(2000). They concluded that the photometric spot effect can nu¬ 
merically reach the fractional spot radius. Then, the feasibility of 
spot detection, studying the spot-induced perturbations in binary 
microlensing events during the caustic crossing and delectabil- 


ity of spot in different wav elengths have been discussed in de¬ 
tails by a number of authors 1 Han et al. 200d : lchang & Han 20021 : 


iHwang & Han 201dl:lHendrv et al. 20021). 

To calculate the magnification factor of a spotted source, we 
assume that the intensity of the spot, i.e. Ip 3 , is almost constant 
over the spot area. Hence, we can re-write the fourth modified 
Stokes parameter S' T (in equation^, by factoring out J/, s as: 


S'r = Sr,* - b I{p 3 ) p\ 

= Sr,*-6 Sr, s , 


/* 1 /*7T 

/ Pdp 

Jo J-1 


dcp A(u) 0(p, (j>) 


(5) 


where we define b = (I(p 3 ) — Ii t3 )/I(p 3 ), I(p 3 ) is the source 
intensity at the geometrical location of the spot center i.e. ( 6 S , 
p B = cos 9 a , Si } + is the Stokes parameter due to the source without 
any spot, given by equation 0 and Sr, s is the Stokes parameter 
of the spot itself considering I(p s ) as its intensity. Integrating is 
done over the source area. In that case, the magnification factor of 
a spotted source is given by: 


Sr,* - b Sr,a 
Sr ,*,0 — b Sr, s ,o 


= A 


1-6/3 

1 - 66 ’ 


( 6 ) 


where the index 0 of the Stokes parameters refers to the ones with¬ 
out lensing effect. A = Sr,*/Sr,*,o is the magnification factor of 
the source without any spot, /3 = Sr, s /Sr,* and 6 = Sr, s ,o/Sr,*,o- 
The magnification factor A' up to the third order in S in terms of un¬ 
perturbed magnification factor A and perturbation terms is given 
by: 

A 1 = A[l + b(5 -p) + 6 2 6(6 -p) + b 3 S 2 {S - /3) + ....]. (7) 


Here, a spot on the surface of source star has two effects: (i) 
decreasing the baseline flux of the source which enhances the 
magnification factor by 66 + 6 2 6 2 + ... and (ii) decreasing the 
magnified light which decreases the magnification factor about 
—6/3 — 6 2 /?6 — .... Before the spot reaches to the caustic curve, 
the magnification factor enhances due to the first effect. When the 
spot light is being magnified due to passing over the caustic line or 
approaching to the lens position, the second effect dominates and 
(4jlhe total magnification factor decreases. When these two effects 
have the same amounts i.e. /3 = 6, the perturbation terms vanish. 
In this time, A' = A = A a where A s = f 1 '’ which means that 

b I,s,0 

the magnification factor of the source is equal to the magnification 
factor of the spot itself. This time does not depend on the relative 
difference in the intensity of source and spot i.e. 6 and just depends 
on the spot area and its location with respect to the lens or caus¬ 
tic line. If the source radius and the lens position (or the caustic 
configuration in binary lensing) in this time are deduced from the 
photometric observations, we can assign the magnification factor 
of the spot itself which helps us to characterize the spot properties. 


2.2 Perturbation on the polarization degree 

Here, we investigate the spot-induced perturbations on polarization 
degree during a microlensing event, some prospects are explained 
in the following. 

(a) Berdyugina (2005) noticed that there are some correlation 
relations between the spot parameters such as its magnetic field, 
temperature contrast and the source photosphere temperature. She 
plotted the data of magnetic field measurement B(kG) and the spot 
temperature contrast A T(K) versus the photosphere temperature 
T*(K) for some active dwarfs and giants and concluded that the 
hotter giants have the weaker magnetic field whereas their spots 
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have the higher temperature contrast. We can fit two useful func¬ 
tions to her data which are given by: 

A T(K) = 1.26 x 10 _4 r* - 0.53T* + 787.40, 

B(kG) = -1.06 x 1CT 3 T* + 7.46. (8) 

According to these relations, for a source star with a constant pho¬ 
tosphere temperature, the cooler stellar spots have the stronger 
magnetic fields and as a result the higher polarization signals. 

On the other hand, the Stokes intensities of the spot, i.e. U iS 
(for i = Q,U, V, I), are proportional to the continuum, local inten¬ 
sity of the spot, i.e. I c , (for more details see Appendix section). But, 
in the microlensing observations from the sources in the Galactic 
bulge, we receive the total flux of the source and its spot. Hence, the 
darker spots have the smaller contributions to the total flux whereas 
they intrinsically have larger polarization signals due to the stronger 
magnetic fields. Here, we investigate which of these two effects, 
i.e. (1) the spot magnetic field and (2) the spot flux contrast, has 
more contribution in the polarimetric perturbations in microlensing 
events. 

For this aim, in Figure |2(a)| we show a polarimetric microlens¬ 
ing event affected by the source spot. We consider a giant star as 
the source with the parameters TfK) = 4314, M*(M 0 ) = 2.6 
and = 24.9. Also we assume that the source has a spot 

with the parameters A T{K) = 831 and B(kG) = 2.9, chosen ac¬ 
cording to the equations {8j. In addition, we consider a cooler spot 
with a stronger magnetic field and the hotter spot with a weaker 
magnetic field. The light curves and polarimetric curves for these 
three cases are shown in left and right panels of this figure. The 
source (black circle) and its spot located over the source surface, 
caustic curve (black curve) and source center trajectory projected in 
the lens plane (red dash-dotted line) are shown with an inset in the 
left-hand panel. The thinner lines in right panel show the intrinsic 
polarization signals of the spotted source. The simple models with¬ 
out spot effect are shown by red solid lines. The photometric and 
polarimetric residuals with respect to the simple models are plot¬ 
ted in bottom panels. Note that the polarimetric residuals are the 
residual in the polarization degree without considering the polar¬ 
ization angles i.e. A P and the absolute value of the residual in the 
polarization vector i.e. | AJ?| whose definition is brought in equa¬ 
tion ©. The parameters used to make this microlensing event can 
be found in Table Q. To calculate the magnification factor of each 
element over the source star, we use the adaptive contouring al¬ 
gorithm JPominik 20071) . According to this figure three spots with 
different magnetic fields and temperature contrasts have almost the 
same intrinsic polarization signals and as a result the same polari¬ 
metric deviations in the polarimetric microlensing. Hence, these 
two effects whereas act reversely have the same strengths. The pho¬ 
tometric signature of the spot depends on the temperature contrast 
and enhances with increasing it. 

The polarization signal of spot maximizes when the spot is 
crossing the caustic curve. If the spot is located on the source edge 
and the source is entering into the caustic from that edge (near the 
spot), the spot-induced perturbation is so higher than when that 
source is coming out from that edge. Because, in the first case the 
light from the other points over the source surface is not magni¬ 
fied. Accordingly, in the first caustic-crossing of Figure |2(a)[ the 
spot-induced perturbation is higher than that in the second one. 

Also, the duration of the polarimetric signal depends on the 
location of spot with respect to the caustic line, the size of spot and 
its host star and it does not depend on the temperature contrast and 
the magnetic field. 

(b) The effect of the projection angle 9 S on the spot-induced 


polarimetric deviation: the amount of the circular Stokes intensity 
in a spectral line due to Zeeman effect is proportional to cos 9 S and 
as a result increases by decreasing 9 S , whereas the linear Stokes in¬ 
tensities in a spectral line are proportional to sin 2 9 S and decrease 
by decreasing 9 S (see equations IA8I and |A5I in Appendix section). 
Hence, the polarization signals due to face-on spots are more cir¬ 
cular and those due to spots with larger projection angles are more 
linear. As mentioned in the Appendix section, the broadband circu¬ 
lar Stokes intensity is so smaller than that in a spectral line by one 
order of magnitude whereas the linear ones have almost the same 
amounts as those in a spectral line. Because, the circular polariza¬ 
tion in a spectral line has an anti-symmetric shape with respect to 
the wavelength. Consequently, stellar spots with the smaller projec¬ 
tion angle 9 S have lower polarized Stokes intensities. 

In Figure [2(b)1 a polarimetric microlensing event of a spotted 
source is represented with three different amounts of the projec¬ 
tion angle 9 S . According to this figure, the intrinsic polarization 
signal due to the spot without lensing for 9 S = 10° (about 0.005 
per cent) is so smaller than that with 9 S = 45° (about 0.026 per 
cent), shown with thin, straight lines. On the other hand, increasing 
the projection angle 9„ decreases the spot area and the total inten¬ 
sity considering limb-darkening effect and as a result increases Sj 
(which decreases the polarization signal). Therefore, by increas¬ 
ing the projection angle, the polarization signal due to stellar spots 
without lensing increases and then decreases. The polarization sig¬ 
nal without lensing due to the spot with 9 S = 80° is equal to 0.018 
per cent, lower than that with 9 S = 45°. 

However, the deviations in the polarimetric microlensing due 
to the spot with 9 3 = 80° is a bit larger than that due to the spot 
with 9 3 = 45° which is due to the finite size effect. The larger pro¬ 
jection angle, the smaller area of the spot. In microlensing events, 
finite size effect always shrinks sharp signals. 

In the photometric curves, increasing the projection angle 9 s 
decreases the spot-induced photometric signal. Because, the spot 
area decreases by increasing the projection angle by cos 9 S which 
decreases the spot-induced perturbation signal. However, the total 
intensity decreases from the center to the limb of the source star due 
to limb-darkening effect (see FigureQJ, but this effect is so smaller 
than the first one. If the area of the spot is fixed by considering the 
limb-darkening effect, the photometric signal due to the spot does 
not change significantly (see e.g. Heyrovsky & Sasselov (2000)). 
Therefore, the spots located at the source center make higher photo¬ 
metric signals as shown in Figure [2(b)1 This points was also noticed 
by Hendry et al. (2002). 

(c) There is a partial degeneracy for indicating the spot size 
and its temperature difference with respect to the host star from the 
photometric deviations. However, these two factors are not com¬ 
pletely degenerate. Because by increasing the spot size the finite 
size effect becomes important. This effect strongly alters the shape 
of the spot-induced signals. Whereas, increasing the spot temper¬ 
ature contrast only increases the perturbation signals and does not 
change the shape of the signals. Hence, these two effect are par¬ 
tially degenerate in the photometric measurements. Although, two 
degenerate spots with the same photometric deviations have dif¬ 
ferent, intrinsic polarization signals, but by considering two dif¬ 
ferent magnetic fields, they can have the same polarization sig¬ 
nal and be polarimetrically degenerate. In Figure |2(c)| we make 
two microlensing events with spotted sources in which the spots 
have the same positions but different sizes, temperature differences 
and magnetic fields. These spots have the same, intrinsic polariza¬ 
tion signals. These two events have the same light and polarimet¬ 
ric curves and are degenerate. However, in the photometric curve 
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Table 1. The table contains the parameters used to make the microlensing events shown in Figures |2(a)||2(b)||2(c)||2(d)| and m respectively. These parameters 
are the source mass M*(M 0 ), the source radius R r (R @), the photosphere temperature of the source TV (K ). the projected source radius in the lens plane, 
normalized to the Einstein radius p+, the spot radius normalized to the source radius r s (R r ), the spot temperature T a (K), the spot magnetic field B{kG), 
the spot projection angles 6° and 0°, the impact parameter of the source center u C m, the angle between the trajectory of source star and the binary axis £°, 
the distance between two lens normalized to the Einstein radius d and the mass ratio of two lens q respectively. The giant stars as sources are chosen from a 
synthetic distribution for giant stars made from Besancon model (Robin et al. 2003). For these figures, we also set the mass of the primary lens AT/ =0.3 Mg, 
the lens and source distance from the observer l)\ = 6.5 kpc and D a = 8.0 kpc and the limb-darkening coefficients cj = 0.7 and ci = 0.04. 


degeneracy is between the spot size and its temperature contrast. 
Therefore, we can not uniquely extract the spot properties from the 
photometric or polarimetric deviations. 

(d) We investigate the finite size effect on the photometric and 
polarimetric signals due to the spot. The spot-induced perturbations 
decrease and lengthen with increasing the source radius. Because, 
by integrating over the source surface, the sharp signals due to the 
points located at the caustic line shrink. In Figure |2(d)| we represent 
a microlensing event with three different amounts of p* and consid¬ 
ering r s = 0.3 R+. The residual for each source size is calculated 
from the spotless source of that particular size. 


2.3 Perturbation on the angle of polarization 

The polarization angle <j> p strongly changes due to the stellar spot. 
Because, this polarization angle is not related to the total Stokes 
parameter, i.e. S' T which is very larger than the polarized Stokes 
parameters. On the other hand, the polarization angle 9 P alters a 
bit. Because, the broad band circular polarization of spots is, one 
order of magnitude, smaller than the linear ones. The spot-induced 
perturbation in the polarization angle causes the polarization vector 
alters as follows: 

| AP| = [P' 2 +P 2 -2P'Pcos6i;cos(</>;-<j!> p )] 1/2 . (9) 

Indeed, we define the polarization vector as P = 
(P cos 6 P cos (p p , P cos dp sin <f> p , Psin 6 V ). Note that in nor¬ 
mal microlensing events without spot effect 9 P = 0. In Figure «[3]l 
we show the curves of the polarization angle (f> p (left panel) and 
polarization degree (right panel) of a microlensing event. 

Measuring the polarization angles helps us to indicate the dif¬ 
ferent components of the polarized intensity which gives some in¬ 
formation about the projection angles 6 a and <j> B . On the other hand, 
measuring the polarization angles in addition to the polarization 
degree can increase detectability of stellar spots through polari¬ 
metric measurements, e.g. in Figure |2(bj| although the maximum 
amount of the residual in the polarization degree A P due to the 
spot reaches to 0.04 per cent, but by considering the variation of 
the polarization angles, the residual in the polarization vector rises 
up to 0.35 per cent. 



Figure 3. A polarimetric microlensing event due a spotted source: The 
curves of polarization angle cf> p and the polarization degree are shown in the 
left and right panels. The simple models without spot effect are shown by 
red solid lines. The polarimetric residuals with respect to the simple models 
are plotted in the bottom of each panel. The parameters used to make this 
microlensing event can be found in Table m. 


3 OBSERVATIONAL PROSPECTS 

In the previous section we studied the spot-induced perturbations 
during caustic-crossing in binary microlensing events. Now, we es¬ 
timate the probability of detecting the source spots in polarimetric 
or photometric observations during caustic crossing. In this regard, 
we perform a Monte Carlo simulation by generating many syn¬ 
thetics caustic-crossing microlensing events with spotted and gi¬ 
ant sources. By considering two useful criteria for polarimetric and 
photometric observations, we survey if the spot-induced perturba¬ 
tions can be discerned in each light curve or polarimetric curve. 
Finally, we estimate the number of events with detectable spot- 
induced signals for given amounts of background stars and obser¬ 
vational time. 

Our criterion for detectability of source spots in the microlens¬ 
ing light curves is ( A 1 — A)/A ^ 2 per cent. For polarimet¬ 
ric measurements, we assume that these observations are done by 
the FOcal Reducer and low dispersion Spectrograph (FORS2) po- 
larimeter at Very Large Telescope (VLT) telescope. This polarime- 
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ter can achieve to the polarimetric precision about 0.1 per cent with 
one ho ur exposure time from a source star brighter than 14.5 mag 
(flngrosso et al. 2015h . To evaluate the detectability of the polari¬ 
metric signals due to the source spots, for each microlensing event 
we evaluate signal to noise ratio (SNR) when the spot is crossing 
the caustic curve, which is given by: 

SNR = J[ 10-°- 4m * A' + n P sFl0- o - 4 ^ fc «]t e:C pl0 0 - 4m *F, (10) 

where, m+(mag) is the apparent magnitude of the source star, A' 
is the magnification factor of the spotted source, Qpsf is the area 
of the Point Spread Function (PSF), p s k v = 22.5 {mag/arcs 2 ) 
is the sky brightness, t exp {s) is the exposure time and m zp = 
28.3 (mag) is the zero point magnitude of FORS2 polarimeter at 
VLT. We assume that there is no blending effect. For detecting the 
spot signal, the exposure time should be less than the time when 
the spot is crossing the caustic curve and we set t exp = tEps/3 
where p g is the spot radius projected in the lens plane and nor¬ 
malized to the Einstein radius. The threshold amount of SNR 
to achieve the polarimetric precision about 0.1 per cent is about 
SNRthr = 34000. Hence, our criteria for detectability of the 
source spot in the polarimetric curves are |AP| >0.1 per cent 
as well as SNR > SNRthr- 

Here, the distribution functions of the binary lenses, source 
and its spot parameters used to generate syntectic binary microlens¬ 
ing events are illustrated. We explained the distribution functions 
for the mass of the lenses, the Galactic coordinates and velocities 
of both sources and lenses, the mass ratio for the binary lenses and 
their distance as well as the sourc e trajectory with respect to th e 
binary axis in our previous works l iSaiadian 2014l : lSaiadian 2013 ). 
and do not repeat them here. 

We consider only giant stars as sources. For indicating the 
mass, surface temperature, absolute magnitude and radius of gi¬ 
ant sou rces, we make a n ensemble of giant stars using Besanqon 
model jRobin et al. 20031) . Then, we choose the source stars from 
that ensemble randomly. We obtain the apparent magnitude of the 
source star (m*) according to the absolute magnitude, distance 
mo dulus and ex tinction. In that case, the more details can be found 
in l lSaiadian 20151) . Also, I consider constant amounts for limb- 
darkening coefficients of source stars as ci = 0.7 and C 2 = 0.04. 

For the spot, we indicate its magnetic field and temperature 
contrast using the equation |8} and adding Gaussian fluctuations. 
There is a correlation relation between the magnetic field of spots 
and the filling factor, i.e. the total fract io nal area co vered by the 
spot, as / = 0.4 B(kG ) — 0.7 l lBerdvugina 20051 ). In the sim¬ 
ulation we consider only one spot over the source surface. For the 
Galactic bulge stars, only the biggest spot can probably be detected. 
We need the distribution function of the radius of the biggest spots 
over source stars. I assume that the ratio of the area of the biggest 
spot to the total area due to all spots over the Sun is the same as that 
ratio due to the spots over other stars. This ratio for the Sun spots is 
about 0.04 dSolanki 20031) . Therefore, I first indicate the total frac¬ 
tional area due to the spots, i.e. /, according to the magnetic field 
by adding a Gaussian fluctuation and then multiply it by 0.04 to 
indicate the fractional area due to the biggest spot. 

In order to obtain the spot position in the lens plane, we need 
two rotation angles of 9 S around y-axis and <p s around z-axis. (p 3 
is uniformly taken in the range of [0,7r]. We choose 6 S by fixing 
the angular surface element. Hence, we uniformly select ip in the 
range of [—1, 1] where 9 S = cos~ lr tp. For the broadband Stokes 
intensities due to the magnetic field, we first calculate them in a 
spectral line and then by considering the factor a (introduced in the 



SNR n 

Figure 4. The scatter plot of the spot-induced polarimetric signals in syn¬ 
thetic microlensing events versus SNR which is normalized to the thresh¬ 
old amount. The blue solid line separates the microlensing events with de¬ 
tectable polarimetric signatures due to the spot. The black stars represent 
the microlensing events with detectable photometric signals due to the spot 
and red circles show the events in which the spot-induced signatures can 
not be discerned through photometric observations. 


Appendix section) estimate the amounts of these Stokes intensities. 
Finally, we investigate with how many probability lensing magni¬ 
fies these polarization signals and makes them be observed. 

Having performing the Monte Carlo simulation, we concluded 
that in about 11.2 per cent of the simulated events the spot-induced 
polarimetric signals are detectable according to our criteria. Note 
that the real efficiency is certainly less than this number. Indeed, 
in observations of microlensing events, 3-5 consecutive data points 
with deviations more than 2-3 ct with respect to the simple model 
are needed to confirm a perturbation effect, where a is the obser¬ 
vational accuracy (see e.g. Dominik et al. 2010). But in this Monte 
Carlo simulation, we do not simulate the data points. Hence, short- 
duration spot-induced signals (e.g. see Figure [2(c)} which are de¬ 
tectable according to our criteria would be difficult to detect in 
real observations due to data spacing and may be confused with 
the noise. On the other hand, observers should fit many synthetic 
light and polarimetric curves with different perturbations to the data 
points. Therefore during observation the true event can be confused 
with other microlensing events which are well fitted to the data 
points e.g. the perturbation due to the stellar spots can be mistaken 
with a transit exo-planet around the source star. 

In Figure 0. we represent the scatter plot of the spot-induced 
polarimetric signals, i.e. |AP| for synthetic microlensing events 
versus SNR normalized to the threshold amount i.e. SNR n . The 
blue solid line separates the microlensing events with detectable 
polarimetric signatures due to the spot from the others. _ 

According to the OGLE-III data l lWvrzvkowski et al. 20141) . 
about /i = 18.4 per cent of the microlensing sources are giant 
stars and about = 5.4 per cent of microlensing events are due 
to the binary lenses with caustic-crossing features. Generally, all 
giant stars are not magnetically active. RS CVn-type binaries, FK 
Coma e-type giants and Lithium rich giants have high magnetic ac¬ 
tivity dKorhonen 20131 ). Recently for some of other giant stars mag¬ 
netic fields have been detected and even mapped (e.g. Strassmeier 
et al. 1990, Auriere et al. 2008, Dorchj2004). Indeed, the magnetic 
activity decreases with the age dSkumanich 19721) . Consequently, I 
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assume that about / 3 = 1 per cent of giant stars have stellar spots 
and magnetic field stronger than 100 G. The probability of detect¬ 
ing a spot over a giant source through polarimetric measurements 
in caustic-crossing of microlensing events can be given by: 


Ts = tfihhT, (IT) 


where t = 0.112 is the polarimetric efficiency for detecting spot- 
induced signals in the simulated events, r is the optical depth for 
the microlensing events towards a given direction. For the direc¬ 
tion of the Galactic Bulge r = 4.48 x 10 -6 JSumi et al. 2006i ). 
hence the optical depth for the polarimetric spot detection is r 3 = 
4.99 x 10 -11 . Finally, we estimate the number of events for Nb g 
background stars being detected during the observational time T 0 b s 
as 


T _ rr T 0 b s Nbg 

s ~ 2 < t E > 


( 12 ) 


where the average Einstein crossing time for the Galactic Bulge 
eve nts d ue to the giant s ources is about < tE >= 24.6 days 
dWvrzvkowski et al. 20141) . By monitoring 150 million objects dur¬ 
ing 10 years towards the Galactic bugle, the number of binary mi¬ 
crolensing events with spotted and giant sources, which have de¬ 
tectable spot-induced polarimetric signals, is 1.75. Finally, we in¬ 
vestigate in how fraction of these simulated events the photometric 
spot-induced signal (A' — A) IA is greater than 2 per cent. In this 
regard, the photometric efficiency for detecting signatures due to 
the spot is 52.6 per cent. In that case the number of events during 
10 years observations from 150 million objects will be 8.19. Note 
that, we consider only giant stars as sources and do not consider 
main-sequence stars. In Figure Q the black stars show these events 
which have detectable photometric signatures due to the spot. 

Consequently, photometry observation is more efficient in de¬ 
tecting stellar spots than the polarimetry observation. However, po- 
larimetry observations during microlensing events give us some in¬ 
formation about the magnetic field of the spots over the sources 
whereas the photometric signals due to the spots do not depend on 
it. On the other hand, the stellar spots change polarization orien¬ 
tation in addition to the polarization degree whereas in the photo¬ 
metric observations of spotted sources only the magnification fac¬ 
tor alters. Variation of the polarization angles give us extra infor¬ 
mation about the different components of the polarized intensity 
and as a result the position of spot on the source surface. Also, the 
polarimetry observations can be complementary to the photometry 
observations for detecting stellar spots with larger amounts of the 
projection angle 9 S . 


stellar spot perturbs polarimetric and photometric curves in binary 
microlensing events during caustic-crossing. 

The stellar spot perturbs the micr olens ing li ght cu rves by mak¬ 
ing positive and negative deviations iHwang & Han 201(1) . If the 
time when the photometric deviation due to spot is zero (between 
positive and negative deviations) is inferred from microlensing 
light curves, we can indicate the magnification factor of the spot 
which helps us to characterize the spot properties except the tem¬ 
perature contrast. 

We figured out some points about the spot-induced perturba¬ 
tions on the polarization degree explained in following: 

(a) The cooler spots have the higher magnetic fields and larger 
polarization signals, but they have the smaller contributions in the 
total flux. We showed that these two properties of spots, tempera¬ 
ture contrast and magnetic field, have the same contributions in the 
polarimetric signal so that increasing both together does not change 
the polarization signal. 

(b) The maximum deviation in the polarimetry curve due to 
the spot happens when the spot is located at the source edge and 
source is entering to the caustic curve from that edge (near the 
spot). 

(c) There is a (partial) degeneracy for indicating the spot size, 
its temperature contrast with respect to the host star and its mag¬ 
netic induction from the deviations in light or polarimetric curves. 

(d) The spot-induced deviations decrease and lengthen by en¬ 
hancing the finite size effect. 

(e) The stellar spots alter the polarization degree as well as 
strongly change its orientation. Variation of the second gives some 
information about the spot projection angles. 

Finally, we estimate the number of the source spots that can be 
observed with this method towards the Galactic bulge. During 10 
years monitoring of 150 million stars towards the Galactic bulge, 
we can specify in the order of 2 stellar spot over giant sources with 
this method whereas about 8 stellar spots over giants can be con¬ 
firmed through photometric observations. Hence, photometry ob¬ 
servation is more efficient in detecting stellar spots than the po¬ 
larimetry observation. However, polarimetry observations during 
microlensing events give us some information about the magnetic 
field of the spots over the sources whereas the photometric signals 
due to the spots do not depend on it. 
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4 CONCLUSIONS 

Gravitational microlensing with finite size effect provides a tech¬ 
nique to probe the anomalies over the source surface specially in 
caustic crossing. One of these anomalies is the stellar spot which 
causes the light intensity of source star does not obey from the 
limb-darkening function. It creates some perturbations in the mi¬ 
crolensing light curves. 

The stellar spots also create a net polarization for the source 
star. Although most spot-induced polarimetric signals of Galactic 
bulge sources are too weak to be detected, but the lensing effect 
can magnify the spot-induced polarization signal and makes it be 
observed. However, lensing effect itself generates a net polarization 
for the source star. Here, we investigated the possibility of detecting 
the source spot through polarimetric observations during caustic¬ 
crossing in microlensing events. In this regard, we studied how the 
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APPENDIX A: THE BROADBAND STOKES INTENSITIES 
DUE TO THE MAGNETIC FIELD OF SPOT 

The aim of this appendix is to estimate the amount of polarized 
Stokes intensities of spots with the magnetic field and take into ac¬ 
count this intensities in the polarimetric microlensing calculations. 
A stellar spot produces a polarization signal through the Zeeman 


effect due to its magnetic field on the spectral lines. The polarized 
radiation can be obtained by solving the Stokes transfer equation: 

— -j— = KI — J , (Al) 

dr c 

where r c is line of sight continuum optical depth, I is the Stokes 
vector, K is the total opacity matrix and is given by: 

K = 1 + vo®, (A2) 

and J is the total emission vector, 


J = ( S c + r/oSi&)eo, 


(A3) 


where 1 is the unit matrix, S c is the source function in the un¬ 
polarized continuum spectra and eo = (1,0,0,0) . We take the 
Minle-Eddington model for the stellar atmosphere where the line 
source function Si = Bo + B\t c and Bo + B\ = I c and I c is 
the continuum intensity close a spectral line. Here, vo is the ratio 
between a spectral line and continuous absorption coe fficien ts. 
is the line absorption matrix which is given by jRees et al. 19891) : 


/ 

VI 

VQ VU 

W \ 


VQ 

Vi Pv 

Pu 


VV 

- pv Vi 

PQ 

V 

Vv 

pU ~ PQ 

m ) 


where 

1 2 1 2 
Vi = -^fesin 9s + -fir + 77 & )(1 + cos 0 a )], 

11 2 

VQ = r-+»?&)] sin 9 S cos 2cf> a , 

11 2 

Vu = -^{Vp --^(Vr + V 6)] sin 0 s sin2<j!> s , 

Vv = IjlVr - Vb) cos 9 S , 

11 2 

PQ = o\pp - o(pr + pb)} Sin 9 a cos2(/)s, 

11 2 

pu = 7j[/°p - o(pr +P6)]sin 9 S sin 2(/> s , 

pv = \{pr - pb) cos 9 S , (A5) 

where 9 a and (f> s represent the position of spot on the surface of 
source. Here the absorption and anomalous p p , r ,b dispersion 

profiles are given by: 


Vp = Vo H(a,v), 

Vb,r = Vo H(a,v±v H ), 
p p = 2vo F(a,v), 

p br = 2vo F(a,v±v H ), (A6) 


in which a is the damping constant, H(a, v ) and F(a,v) are Voigt 
and Faradey-Voigt profiles, v = (A — Ao)/AAd and vh = 
AAz/AAd are the wavelength separation from the spectral line 
center (i.e. Ao) and the wavelength shift due to the Zeeman ef¬ 
fect. Both parameters are normalized to the Doppler broadening 
(i.e. AA_d). The second one is given by: 

AAz = eA . 2i?ge// , (A7) 

An m c 

in which B is the value of the magnetic field, g e f / is the effective 
Lande factor and e, m and c have their usual meanings. In this 
work we used the averaged values of these parameters: a = 0.1, 
AAd = 300 mA, vo = 1 0, A = 60 00A, g e ff = 1.5 and Bo = 
B\ = 0.5 I c jHuovelin & Saar 199 ill . 
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By assuming some limitations for the stellar at mosphere , 
there are analytical solutio ns for the Stokes in ten sities dlJnn o 19561 : 

lAuer et al. 19771 : lLandolfi & Landi Deerinnocenti 1982^) . Bv con- 

sidering the following limitations: (i) the magnetic field is con¬ 
stant in amount and direction over the spot, (ii) linear de¬ 
pendence of the source function with optical depth (i.e. Bo 
and B i) and (iii) constant ratio of the line and continu- 
ous absor bt ion coefficients, the Stoke s intensities are given by 
( lLandolfi & Landi Desil'Innocenti 1982h : 

I\ = Bo + fiBiA x [(l + ?7-r)((l + Vi) 2 + Pq + Pu + Pv)], 

Q\ = -p BiA _1 [(l + rn) 2 riQ + (1 + r/i)(wPu ~ WPv) + 

+ Pq(vqPQ + WPu + WPv)], 

U\ = -p B 1 A~ 1 [(l+rp) 2 w + (l + Vi)(VQPv ~ WPq) + 

+ Pu(VQPQ + WPu + WPv)], 

V\ = -p BiA _1 [(l + pi) 2 w + 

+ pviVQPQ + Vupu + WPv)], (A8) 

where the index of A refers to the Stokes intensities in a spectral 
line and 

A = (1 + »?/)"[(l + r)i)~ — tiq — rpj — r/y + p 2 Q + 

+ Pu + Pv] - VlQPQ + Vupu + WPvf- (A9) 

The microlensing observations are done in optical band and 
we should estimate the contribution of the spectral lines in the 
broadband photometric mode observation. In this regard, we as¬ 
sume that the broad band polarization can be approximated by us¬ 
ing the polarization in an average line profile. Then, we scale the 
resulted Stokes intensity with a factor to obtain the Stokes intensi¬ 
ties in a given pass band. This method was first used to calculate 
the broad band linear polarization in cool stars, (e.g. Calamai et al. 
(1975), Huovelin & Saar (1991), Saar & Huovelin (1993)). Hence, 
we compare the amounts of the spectral Stokes intensities with the 
amounts of broadband for fixed amounts of the magnetic field and 
the projection angle. Indeed, we assume the similar dependence of 
the broadband Stokes intensities to angles of 6 S and <j) 3 while keep¬ 
ing same magnetic field. 

We represent the broadband Stokes intensities due to the spot 
magnetic field (equation[4]i with Ii t3 = i\a (for i = I,Q, U, V), in 
which a is the ratio of the broadband Stokes intensities to the those 
in a spectral line. 

Considering Milne-Eddington model for the source atmo¬ 
sphere, for the case of B = 3.2 kG and 0 = 0°, the maximum 
amo unt of t he broadband circular polarization reaches to 0.1 per 
cent dStift 1996ll . while the circular polarization in spectral line of 
630lA is about 3 — 4 per centQ. Therefore, the broadband circular 
polarization is one order of magnitude less than the circular polar¬ 
ization in a spectral line. Indeed, the Stokes V profile has approxi¬ 
mately an anti-symmetric shape versus v in a spectral line. Hence, 
the net circul ar polarizatio n over a broad band would be negligi¬ 
ble dBorrero & Ichimoto 201 ll ). For B = 4 kG and 6 = 60°, the 
broadband lin ear a nd circular polarizations are about 0.8 and 0.02 
per cent dStift 19971) . For the spectral line of 6301A they are almost 
1 and 0.7 per cent. Comparing these numbers, the broadband linear 
polarization is almost equals to the linear polarization in the spec¬ 
tral line. Therefore, we set a = 1 for linear Stokes intensities and 
a = 0.1 for the circular Stokes intensity. 


1 http://www.iac.es/proyecto/magnetism/pages/codes/milne-eddington- 
simulator.php 













